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Abstract 

We study and compare the persistence of bipartite and multipartite entanglement in one and 
two-dimensional spin XY model in an external transverse magnetic field under the effect of thermal 
excitations. Wc compare the threshold temperature at which the entanglement vanishes in both 
cases. We use the concurrence as a measure of the bipartite entanglement and the geometric 
measure to evaluate the multipartite entanglement of the system. We have found that for the 
anisotropic and partially anisotropic systems the nearest neighbor bipartite entanglement vanishes 
asymptotically at much higher magnetic field compared to both the next to nearest neighbor 
bipartite entanglement and the multipartite entanglement which asymptotically coincide. Also the 
same behavior was observed for the threshold temperatures where the nearest neighbor bipartite 
one is much higher than both of the next to nearest neighbor bipartite and multipartite where the 
latter two coincide asymptotically and the three of them increase monotonically with the magnetic 
field strength. Thus as the temperature increases to certain value, the threshold, the multipartite 
entanglement and the bipartite entanglement of the far parts of the system may vanish while the 
nearest neighbor bipartite entanglement may sustain up to much higher temperature. For the 
isotropic system, all types of entanglement and threshold temperatures vanish at the same exact 
small value of the magnetic field. We emphasis the major role played by both the properties of the 
ground state of the system and the energy gap as well. Furthermore, we found that the quantum 
effects in the system can be maintained at high temperatures, where we have observed that the 
different types of entanglements in the lattice sustain at high temperatures if we apply sufficiently 
high magnetic fields. 

PACS numbers: 03.67.Mn, 03.65.Ud, 75.10.Jm 
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I. INTRODUCTION 



The state of a classical composite system is described in the phase space as a product of 
its individual constituents separate states but the state of a composite quantum system is ex- 
pressed in the Hilbert space as a superposition of tensor products of its individual subsystems 
states. Therefore the state of a quantum composite system is not necessarily expressible as 
a product of the individual quantum subsystems states. This peculiar property of quantum 
systems is called Entanglement, which has no classical analog Recently the interest in 
studying quantum entanglement was sparked by the development in the fields of quantum 
information and quantum computing which was initiated in the eighties by the pioneering 
work of Benioff, Bennett, Deutsch, Feynman and Landauer 2|-[8|. Although there is still no 
complete theory that can quantify entanglement of a general multipartite system in pure 
or mixed state, there are few cases where we have successful entanglement measures. Most 
importantly, bipartite system in a pure state and mixed state of two spin 1/2, also for pure 
and mixed multipartite systems using geometric measures such as geometric entanglement 



and relative entanglement [9| 



14| . Quantum information processing and quantum computa- 



;ions can only be performed in a many body system with very sophisticated arrangements 



15(1 . The building unit, smallest for storing information in such a system (qubit), has to 



be a well defined two state quantum entity that can be easily addressed, manipulated and 
readout. The basic idea is to define certain quantum degree of freedom to serve as a qubit, 
such as the charge, orbital or spin angular momentum. The next step is to define a control- 
lable mechanism to form coupling between two individual qubits in such a way to produce a 
fundamental quantum computing gate. In the meantime, we have to be able to coherently 
manipulate such an entangled state to provide an efficient computational process. Natural 
systems of interest have strong interaction with its environment, which makes it a signifi- 
cantly challenging mission to achieve the high coherence control required to manipulate the 
system, where decoherence is considered as one of the main obstacles toward realizing an 



effective quantum computing system 16l-ll9[|. Particularly, practical systems are required 



to function at non-zero finite temperatures, which means that the system will be exposed 
to thermal excitations and therefore its mixed thermal states should be fully studied and 
understood. Evaluating the density matrix of mixed thermal states of many body systems 
is very hard task due to the large size of the Hilbert space of the system in that case. Nev- 
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ertheless, there has been an approach to provide a transition temperature below which the 

multipartite entanglement is guaranteed in such systems based only on information about 

I I 

the ground state of the system and its partition function [20]. Also recently the robustness 
of thermal multipartite entanglement in one-dimensional XY system was investigated and 
the corresponding transition temperature was estimated Quantum phase transitions 

in many body systems are accompanied by a significant change in the quantum correlations 
within the system, which led to a great interest in investigating the behavior of quantum 



entanglement close to the critical points of transitions |22|-|25| . This raised the question of 
the multipartite versus bipartite entanglement and whether they have to coexist and which 
one of them is more capable of describing the critical behavior of many body systems. It was 
demonstrated that there may exist a mixed state which contains no bipartite entanglement 
but possesses multipartite entanglement [2^ . 

In this paper, we consider two different systems of finite number of spins, each in presence 
of an external transverse magnetic field in contact with a heat bath at temperature T. The 
two systems are one dimensional XY spin chain and two-dimensional XY triangular spin 
system, where the number of spins in each system is 7 and the nearest neighbor spins are 
coupled through an exchange interaction J. We investigate and compare the bipartite and 
the multipartite entanglement of both systems under the effect of external magnetic field and 
thermal excitations and different degrees of anisotropy. We show that, for both systems, the 
nearest neighbor bipartite entanglement is more resistant to thermal excitations compared 
to the next to nearest bipartite and the multipartite entanglements where the latter two 
coincide asymptotically as the magnetic field strength increases. Also we demonstrate that 
the threshold temperature, at which the entanglement vanishes, is much higher for the 
nearest neighbor entanglement compared to both of the next to nearest neighbor bipartite 
and multipartite entanglements. We also examined the persistence of quantum effects at high 
temperatures by observing the entanglement behavior and show that we may maintain non 
zero entanglement at considerably high temperatures by applying strong enough magnetic 
fields. 

This paper is organized as follows. In the next section we present our model. In sec. 
Ill we focus on the two-dimensional spin system and evaluate the bipartite entanglement 
and the thermal energy of the system. The multipartite entanglement and the threshold 
temperatures for all type of entanglements for the two-dimensional system are evaluated in 
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sec. IV. In sec. V we calculate and compare the bipartite and multipartite entanglements 
and the corresponding threshold temperatures in the one-dimensional system. We conclude 
in sec VI. 

II. THE MODEL 

We consider two different systems, two (and one) dimensional spin-1/2 XY model subject 
to an external magnetic field h. The number of spin in the system is 7, which are localized on 
a two-dimensional triangular lattice (one dimensional chain) with an exchange interaction 
between the spins J as shown in fig. [TK (b). The Hamiltonian of the system is given by 




FIG. 1: (a) Two-dimensional triangular spin lattice; (b) One-dimensional spin chain. 

neighbors sites on the lattice and Jij = J for all sites. For this model it is convenient to 
study a dimensionless Hamiltonian where we set J = 1 and define a dimensionless parameter 
A = h/J. The Hilbert space of this spin systems is huge with 2^ dimensions, nevertheless 
it can be exactly diagonalized using the standard computational techniques, yielding the 
system energy eigenvalues {Ei} and eigenfunctions {ipi}. At absolute zero temperature, the 
system lies in its ground state {ipo) which is usually entangled with an amount that varies 
based on the values of the different system parameters. The system is described by the 
density matrix defined in terms of the pure ground state wavefunction \tpo) as 

P=|V^o)(^o|- (2) 
5 



Now when the spin system is set into contact with a heat bath at an absolute temperature 
T, the system moves from its initial pure state described by eq. ([2]) to a mixed thermal 
state, which is a mixture of the ground state and a number A'^e of excited states, represented 
by 

PT = ^{e-^^lV^o)(^o| +$^e-'^^«|^.)(z^.|} , (3) 

1=1 

where /3 = 1/kT, k is Boltzmann constant and Z is the system partition function. The 
number of excited states involved depends on temperature, where more states are added 
as the temperature is raised. This mixing of excited states with the ground state act as a 
destructive noise that reduces the amount of entanglement contained in the system. When 
the temperature reaches certain value, which varies based on the system characteristics 
and parameters values, the amount of noise created by the excites states due to thermal 
fluctuations is sufficient to turn the system into a disentangled state. This temperature 
is known as the threshold temperature, denoted by Tth , where below it the system is 



guaranteed to be entangled 



20|. 



III. BIPARTITE ENTANGLEMENT IN TWO-DIMENSIONAL SPIN SYSTEM 

To study the bipartite entanglement in the system, we confine our interest to the entan- 
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. All the information about the 



glement between only two spins, at any sites i and j 
considered two sites i and j is contained in the reduced density matrix pij which can be 
obtained from the entire system density matrix by integrating out all the spins states except 
i and j. We adopt the entanglement of formation, as a well known measure of entanglement 
where Wootters [29I] has shown that, for a pair of binary qubits, the concurrence C, which 
goes from to 1, can be used to quantify entanglement. The concurrence between two sites 
i and j is defined as 

C{pij) = max{0, ei - €2 - €3 - €4}, (4) 



where the e^'s are the eigenvalues of the Hermitian matrix R = a/ \fpp\fp with p = (a^ ® 
ay)p*{ay ® cr^) and o"^ is the Pauli matrix of the spin in y direction. For a pair of qubits the 
entanglement can be written as, 

E(p,,) = e(C(p,,)), (5) 
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(b) 






FIG. 2: (Color online) The concurrences C(l,2), C(l,4) and C(l,7) of the two-dimensional Ising system 
(7 = 1) versus A and fcr(in units of J). 



where e is a function of the "concurrence" C 

.(c) . . 

where h is the binary entropy function 

h{x) = —X log2(x) — (1 — s) log2(l — x). 



(6) 



(7) 



The concurrence can be used as a measure of entanglement as it coincides with the entan- 
glement function in behaviour and range. In our calculations we use the concurrence C as 
a measure of the bipartite entanglement. Using the mixed density matrix px defined in 
([2]), one can evaluate the bipartite entanglement between any pair of spins in the system. 
In this section we focus on studying the bipartite entanglement only in the two-dimensional 
triangular XY spin system sketched in fig. W^a). In fig. [2] we have explored the behaviour 
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FIG. 3: (Color online) The contour plot of the concurrences C(l,2), C(l,4) and C(l,7) of the two- 
dimensional Ising system (7 — 1) versus A and fcT(in units of J). 

of the concurrences of the nearest neighbors C(l, 2); C(l,4) and the next-to- nearest neigh- 
bors C(l,7) versus A and the temperature KT for the anisotropic Ising system (7 = 1). 
As can be noticed, the nearest neighbor bipartite entanglements between two border sites 
C(l,2) and between a border site and the central one C(l,4) are strongest for very small 
magnetic field but very fragile away form the zero temperature. On the other hand, as the 
magnetic field is increased the entanglement maintains a small value which is more resistant 
to higher temperatures. Interestingly, the threshold temperature Tfh at which the entangle- 
ment vanishes increases monotonically as the magnetic field increases. The next to nearest 
neighbor entanglement C(l,7) sustains only for very small values of magnetic field and in 
the vicinity of the zero temperature and its value is much smaller than the nearest neighbor 
entanglements. In order to investigate the thermal robustness of the entanglement state 
we show the contour plot of the concurrences C(l,2), C(l,4) and C(l,7) at much higher 
temperatures in fig. [31 As one can see, we can reach, for C(l,2) and C(l,4), a threshold 
temperature higher than kT = 40 by applying a magnetic filed as high as h = 150 with an 
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FIG. 4: (Color online) The concurrences C(l,2), C(l,4) and C(l,7) of the two-dimensional partially 
anisotropic system (7 = 0.5) versus A and A;T(in units of J). 

entanglement value of the order 10^^. The next-to-nearest neighbor concurrence C(l,7) is 
very fragile to temperature regardless of the strength of the applied magnetic field as shown 
in fig. El^c). The partially anisotropic system with 7 = 0.5 was found to exhibit a close 
behavior to the 7 = 1 case as depicted in figs. H] and |5l The peak of the concurrences, at 
small magnetic field and in the neighborhood of zero temperature, is not single and this is 
due to the profile of the system energy gap for 7 = 0.5 as will be discussed in more details 
below. It is clear that the value of the threshold temperature corresponding to the different 
magnetic filed values is smaller compared to that of the 7 = 1 case as can be concluded from 
figs. [3] and [5l Interestingly, the completely isotropic spin system with 7 = behaves in a 
completely different way compared to 7 = and 0.5 as shown in figs. |6]and [71 As can be 
noticed from the figures, not only the next-to- nearest neighbor concurrence C(l, 7) but also 
the nearest neighbor concurrences C(l,2) and C(l,4) vanish at very small temperatures. 
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FIG. 5: (Color online) The contour plot of the concurrences C(l,2), C(l,4) and C(l,7) of the two- 
dimensional partially anisotropic system (7 = 0.5) versus A and fcr(in units of J). 

about few kT. Clearly the thermal fluctuations is very devastating to the isotropic system 
where the bipartite entanglements over the whole lattice vanishes at very small temperature. 
The peak of the concurrence C(l, 2) is higher than that of C(l, 4) but C(l, 7) is much lower 
than both. In fact this behavior of the isotropic system (with 7 = 0) should not be very 
surprising as this system belongs to a different universality class from the one corresponding 
to < 7 < 1, where its ground state is already separable for A > 1.85. In order to investigate 
the robustness of entanglement at much higher temperatures, we depict the contour of the 
concurrences C(l,2) and C(l,4) in the anisotropic system at very high magnetic field in 
fig. m which confirms the survival of entanglement, though very low in the order of 10"'', at 
as high temperature as kT = 130. It is of great interest to examine the relationship between 
the robustness of thermal entanglement and the corresponding thermal energy gap of the 
system. By the thermal energy gap, AEfh, we mean the difference between the mean (en- 
semble average) energy of the system at temperature T and the system ground state energy, 
i.e. AEth = {E) — Eq. In fig. [9]^a) we plot the energy gap of the anisotropic system versus A 
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FIG. 6: (Color online) The concurrences C(l, 2), C(l, 4) and C(l, 7) of the two-dimensional isotropic system 
(7 = 0) versus A and fcr(in units of J). 

and temperature kT, and in figs. IHl^b), (c) and (d) we present the contour plots of energy gap 
for the systems with 7 = 1, 0.5 and respectively. Remarkably, one can see a strong corre- 
spondence between the strength and survival of entanglement, particularly for 7 = 0.5 and 
1 and the value of the energy gap when comparing figs. |3l [5] and [H The energy gap is either 
zero (the white regions of the contour plot) or quit small in the domains of non- vanishing 
entanglement. Two important notices are due here, first, the energy gap profile, presented 
in fig. ini looks almost the same for the three anisotropic parameter values, nevertheless the 
entanglement in the isotropic system (7 = 0) in contrary to the other two cases vanishes at 
very low temperature regardless of the energy gap value and this is due to the fact that this 
system ground state, as we mentioned before, is disentangled for A ~ 1.85 and higher. The 
other important point is that the thermal energy gap at the different anisotropic values looks 
asymptotically (at high magnetic field) the same, but one can see clear differences between 
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FIG. 7: (Color online) The contour plot the concurrences C(l, 2), C(l, 4) and C(l, 7) of the two-duTiensional 
isotropic system (7 — 0.) versus A and fcr(in units of J). 

them at small magnetic field values as shown in fig. [TOl This explains the different profiles 
of the concurrence peaks, at small values of the magnetic field, as the degree of anisotropy 
changes as demonstrated in figs. [3l [5] and [71 

IV. ROBUSTNESS OF THERMAL ENTANGLEMENT IN TWO-DIMENSIONAL 
SPIN SYSTEM 



In order to study the multipartite entanglement of the entire lattice, a distance- 



measure of entanglement, namely the global robustness of entanglement R{p) 20|, l30[ is 



ike 



commonly used, which is defined for a general state p as the minimum amount of noise t 
needed to destroy the entanglement content of p and is given by 



R{p) := min t 



(8) 



12 




300 350 400 450 500 300 350 400 450 500 

I. X 



FIG. 8: (Color online) The contour plot of the concurrences C(l,2), C(l,4) and C(l,7) of the two- 
dimensional Ising system (7 = 1) versus A and fcT(in units of J) for a range of A from 300 to 500. The color 
map is as shown in subfigure (a). 
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FIG. 9: (Color online) (a) The thermal energy gap (in units of J) of the two-dimensional Ising system 
(7 = 1) versus A and fcT(in units of J); (b) the contour plot of thermal energy gap of the two-dimensional 
spin system versus A and fcT(in units of J) for 7 = 1; (c) 7 = 0.5 and (d) 7 = 0. 
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(c) 




FIG. 10: (Color online) The contour plot of thermal energy gap (in units of J) of the two-dimensional spin 
system versus A and fcr(in units of J) for (a) 7 = 1; (b) 7 = 0.5 and (c) 7 = 0. The color map for all 
subfigures is as shown in subfigure (a) 



where u is the state when added to p convert it to a separable state such that 

1 



1 + t 



(9) 



where S is the set of all separable states. The resultant separable state (p can be regarded as 
a mixture of two states p and u with relative populations 1/(1 + 1) and t/{l + t) respectively. 
This general approach can be applied, in particular, to a system in contact with a heat 
reservoir to determine the threshold temperature, T^h, below it the system is guaranteed to 



be entangled 



20|. Therefore, for instance in the spin system, if the state p is identified as the 



ground state ipo and the rest of the states {ipi}, which get mixed with ipQ as the temperature 
is raised, as to, then the population of the state p is given by 1/(1 + t) = e~^°^'''^ /Z. As a 
result, the condition for the system to be guaranteed entanglement at a temperature T will 
read 



^-Eo/kT 



> 



(10) 
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where R{ipo) is the global robustness of the ground state ipo, which has an energy eigenvalue 
Eq. To obtain the threshold temperature one has to turn the inequality in Eq. dlOp into an 
equality and we get 

p-Eo/kTth 1 

= = . (11) 

To determine the threshold temperature, Tth-, one has to evaluate the robustness of en- 
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FIG. 11: (Color online) The concurrences C(l,2), C(l,4), C(l,7) and geometric entanglement of the two- 
dimensional spin system virus A for 7 = 0, 0.5 and 1 at zero temperature. The legends are as shown in 
subfigure (a). 

tanglement of the ground state R{ipo) which is very difficult task, where the entire system 
Hilbert state has to be searched for the noise mixed state u. However, it was found that 
a lower bound for the robustness of entanglement can be obtained j20| by evaluating the 
geometric entanglement G{iI)q) instead which is easier to evaluate, where in general for 
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FIG. 12: (Color online) The threshold temperatures (in units of J) corresponding to the concurrences C(l,2), 
C(l,4), C(l,7) and geometric entanglement of the two-dimensional spin system virus A for 7 = 0, 0.5 and 1 
at zero temperature. The legends are as shown in subfigure (a). 



any pure state 



1 



1 + R^i)) - 2 ' 

which would enable us calculating a lower bound for the threshold temperature, where 
below it the system is guaranteed to be entangled. In all the figures we simply denote this 
temperature as Tth- 

The geometric measure of multipartite entanglement utilizes the geometric properties of 
the Hilbert space to find the distance (or angle) between a pure state, ^Z', representing the 
system and the closet pure separable state, to it, i.e. || IV') — |0) ||- The square sin of the 
angle between the two states ^/^ and represents a good measure of the global geometric 
entanglement, where the smallest value of the square sin specifies the closet separable state 
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FIG. 13: (Color online) The energy gap (in units of J) of the two-dimensional spin system virus A for 
7 = 0, 0.5 and 1 at zero temperature. 
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FIG. 14: (Color online) (a) The geometric entanglement of the two-dimensional Ising system (7 — 1) and 
(b) the threshold temperature (in units of J) of the two-dimensional isotropic system (7 — 0), for two 
different numbers of the set of parameters {Pi)-, n = 2 x 10^ and 4 x 10^. 
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to the pure state ijj and is defined by 

G{^):=l-[me^x^\\{m\\]\ (13) 

where represents cosine the angle between the two states ip and 0. By evaluating 

the the geometric entanglement (GE), and using Eqs. (|TT|) and f|T2l) we can find the lower 
limit of the threshold temperature Tth- In order to find the closet separable state to the state 
ip, utilizing the reality of the wavefunction and the symmetries of the system, we assume 
that the arbitrary separable state is a product of the single spin states of the 7 spins which 
takes the form 

\<p) = llm + ^i'Pm (14) 

i=l 

The set of parameters {Pi,i = 1,2, ■■ - 7} has to be varied over its entire range to cover 
the whole system Hilbert space searching for the closet distance between ip and (p where 
< Pj < 1. Searching the entire Hilbert space is a computationally difficult task and 
therefore we have tested around 4 x 10^ different distinct set of values for the P's parameters, 
uniformally distributed over the system Hilbert space, for the calculations of the geometric 
entanglement. 

In this section we investigate only the two dimensional spin system. In fig. [TT], we com- 
pare the bipartite entanglements C(l,2), C(l,4) and C(l,7) with the multipartite geomet- 
ric entanglement GE versus the parameter A at zero temperature for different degrees of 
anisotropy. As can be noticed in fig. [TTT a). the nearest neighbour bipartite entanglements 
C(l,2) and C(l,4) of the Ising system decay as A increases but vanishes at large values of 
the parameter compared to both the next to nearest bipartite entanglement C(l, 7) and the 
multipartite entanglement GE where they reach zero magnitude at much lower values of 
lambda. This behavior of the bipartite and multipartite entanglement is repeated again in 
the partially anisotropic system as shown in fig. [TlTb) where the asymptotic behavior of the 
nearest neighbor entanglements coincide and that of the next to nearest and multipartite 
entanglements coincide but at smaller values of A compared to the Ising case. The behavior 
of the entanglement in the isotropic system is depicted in fig. [TlTc) where all types of entan- 
glement vanish at A ~ 1.85 and as we mentioned before this is due to the fact that the ground 
state of the system is separable at this value and higher. In fig. [12] we compare the threshold 
temperature of the bipartite entanglements, which is the temperature at which the bipartite 
entanglement vanish as was demonstrated in figs. [3l [5] and [71 to the threshold temperature 
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of the multipartite geometric entanglement, as defined before versus the parameter A. In 
the Ising model, explored in fig. [T2r a). the threshold temperatures of the nearest neighbour 
bipartite entanglements C(l, 2) and 4) coincide and increase monotonically as the mag- 
netic field increases. On the other hand, Tfh for the next to nearest neighbor entanglement 
C(l,7) coincides with that of the geometric entanglement at small values of the magnetic 
field where it increases monotonically but suddenly drops to zero around A = 6, while Tth for 
the geometric entanglement maintains its monotonic behavior but is much smaller than Tth 
for C(l, 2) and 4). In fig. [T2r b). the threshold temperatures of the partially anisotropic 
system, 7 = 0.5, behave in a similar way to the isotropic case where the temperatures for 
the nearest neighbor bipartite entanglements coincide but what is even more interesting is 
that the threshold temperatures for the next to nearest neighbor bipartite and geometric 
entanglement coincide too. This means that at some threshold temperature the multipartite 
entanglement within the lattice vanishes along with the bipartite entanglement between the 
far spins (such as the next to nearest case 7)). Nevertheless the bipartite entanglement 
between the nearest neighbor spins, such as C(l, 2) and C(l, 4) sustains at the same temper- 
ature but may vanish at higher temperature. A closer look at the behavior of the threshold 
temperatures of the partially anisotropic system at small values of the magnetic field is 
given in fig. [TW c). One can see a sharp changes in the threshold temperatures specially for 
the next to nearest bipartite entanglement and the geometric entanglement, which can be 
explained in terms of the energy gap of the system as will discussed next. The threshold 
temperatures of the completely isotropic system, 7 = 0, is explored in fig. [T2r d). where 
again the threshold temperatures of the nearest neighbor entanglements C(l, 2) and C(l, 4) 
are very close and maintain an almost constant value before suddenly dropping to zero at 
A ~ 1.85. The threshold temperature of the next to nearest neighbor bipartite entangle- 
ment, which considerably lower than that of the nearest neighbors, increases linearly before 
suddenly dropping to zero also at the same value A ~ 1.85. The threshold temperature of 
the multipartite entanglement exhibits an oscillatory behavior with an average value within 
that of the next to nearest bipartite value. Figure [13] explains the sharp changes in the 
threshold temperatures at the small range of values of the magnetic field. As can be noticed 
in fig. [TST a). (b) and (c), the energy gap of the Ising system increases monotonically over 
the entire A range with no sharp changes, which explains the smooth monotonic increase in 
the threshold temperatures shown in fig. IT^ a). In contrary, the energy gaps in the partially 
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anisotropic and isotropic systems exhibit oscillating and sharp oscillating changes respec- 
tively at the small values of the parameter A < 1.85 before coinciding with the anisotropic 
curve and increasing monotonically as shown in fig. [T31 Interestingly, by comparing the 
behavior of the threshold temperatures, particularly of the multipartite entanglement, in 
fig. [12] to that of the energy gaps in fig. [T^ one can notice the strict correspondence between 
them; the minima (and maxima) in the threshold temperatures coincide with that of the 
energy gaps on the parameter A scale and when the energy gap increases monotonically, 
the threshold temperature follows that behavior too. The impact of the energy gap on the 
multipartite entanglement is stronger compared to the bipartite entanglement due to the 
fact that the energy gap is calculated for the entire (multipartite) system. The effect of the 
number of distinctive set of parameters = 1,2, ■■ -7} on the accuracy of the results 

is examined in fig. [H] where two different numbers, 2 x 10^ and 4 x 10^ are compared in 
plotting the multipartite entanglement for 7 = 1 and the threshold temperature for 7 = 0.5 
, which shows a very close coincidence. 

V. ENTANGLEMENTS AND THRESHOLD TEMPERATURES IN ONE DIMEN- 
SIONAL SPIN SYSTEM 

Now let us consider a one dimensional XY spin chain consisting of 7 spins, as sketched in 
fig. IH^b), which is described by the same Hamiltonian Eq. ([T]) where in this case the exchange 
interaction Jij exists only between each spin and its two nearest neighbor spins on the chain. 
The system shows a close behavior to what we have seen in the two-dimensional case. In 
fig- [TST a) we compare the bipartite entanglements to the multipartite entanglement for the 
Ising system, where as can be seen the nearest neighbor entanglements C(l, 2) and C(2, 3) are 
very close where they start with zero magnitude at A = and increase as A increases reaching 
a maximum value around A = 1 and then decays to zero at large A. The next to nearest 
neighbor entanglement C(l, 3) exhibits a similar behavior but with much smaller magnitude 
and asymptotically coincides with the multipartite entanglement where both reach zero 
magnitude at smaller A value compared to that of C(l, 2) and C(2, 3). The behavior of the 
partially anisotropic system is very close to that of the Ising system as shown in figs. [TST b) 
and (c) except the quasi-oscillatory behavior of the bipartite entanglement at values of A < 1 
but again the nearest neighbor entanglements coincide asymptotically as well as both the 
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FIG. 15: (Color online) The concurrences C(l,2), C(l,4), C(l,7) and geometric entanglement of the one- 
dimensional spin system virus A for 7 = 0, 0.5 and 1 at zero temperature. The legends are as shown in 
subfigure (a). 

next to nearest neighbor and the multipartite entanglements. The entanglements of the 
isotropic system, similar to the two-dimensional case shows a step-like behavior and vanish 
at the same point, which is A ~ 0.92 in the current case as depicted figs- ITST d). The threshold 
temperature of the different types of entanglements in the one dimensional chain is explored 
in fig. [161 Once more the behavior of the threshold temperatures of the nearest neighbor 
entanglements C(l,2) and C(2,3) are close at the different degrees of anisotropy of the 
system and the case is the same for the next to nearest neighbor bibartite and multipartite 
entanglements and as can be noticed the isotropic system has zero threshold temperature 
at A ~ 0.92. The behavior of both the entanglements and threshold temperatures in the 
one- dimensional spin chain can be explained in terms of the variation in the system energy 
gap at zero temperature which is depicted in fig. [T71 Similar to the two-dimensional case. 
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FIG. 16: (Color online) The threshold temperatures (in units of J) corresponding to the concurrences 
C(l,2), C(l,4), C(l,7) and geometric entanglement of the one-dimensional spin system virus A for 7 = 0, 0.5 
and 1 at zero temperature. The legends are as shown in subfigure (a). 

one can see the strict correspondence between the variations in entanglements and threshold 
temperatures and the variations in the energy gap for all degrees of anisotropy of the system 
and at all A values. 

An estimate of the experimental values of the threshold temperatures for the typical spin 
systems of interest are due here. As the values of the threshold temperatures and energy 
gaps are all expressed in units of the exchange interaction constant J, which varies for the 
considered spin systems over a range of the order of fieV to meV 3ll-l34| , the corresponding 
range of the threshold temperature is 1.16 x 10"^ K to 11.6 K. Also the typical value of the 
magnetic field h, which is also is expressed in units of J, can be evaluated here which goes 
over the range 1.7 x 10^^ T to 17 T. Our results demonstrate that bigger energy gap would 
lead to higher threshold temperature but needs stronger magnetic field too. Using these 
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FIG. 17: (Color online) The energy gap (in units of J) of the one-dimensional spin system virus A for 
7 = 0, 0.5 and 1 at zero temperature. 

results, one can come up with important estimates, where as can be concluded from fig. [121 
the two-dimensional Ising system can reach a bipartite threshold temperature as high as 100 
K which needs a magnetic field that is as high as 300 T but the corresponding multipartite 
threshold temperature would be only about 50 K. The isotropic system is entangled up to 
a magnetic field of about 30 T where the maximum bipartite threshold temperature would 
be about 15 K and the maximum reachable multipartite threshold temperature is 7.5 K. 
For the same applied magnetic field, the threshold temperatures of the one dimensional spin 
chain would be slightly smaller than the corresponding ones in the two dimensional system 
as can be concluded from fig. [161 
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VI. CONCLUSION AND FUTURE DIRECTIONS 



We have investigated the robustness of bipartite and multipartite entanglement in one 
and two-dimensional XY spin lattices in an external magnetic field h against thermal excita- 
tions. The spins are coupled to each other through nearest neighbor exchange interaction J. 
The number of spins in the lattice is 7, which are coupled to a heat bath at temperature T. 
We have compared the bipartite entanglement to the multipartite entanglement versus the 
external applied magnetic field and temperature. Also we compared the threshold tempera- 
ture at which the entanglement vanishes in both cases. We use the concurrence as a measure 
of the bipartite entanglement and the geometric measure to evaluate the multipartite entan- 
glement of the system. We have observed that for the anisotropic and partially anisotropic 
spin systems the nearest neighbor bipartite entanglement vanishes asymptotically at much 
higher magnetic field compared to the next to nearest neighbor bipartite entanglement and 
the multipartite entanglement where the latter two asymptotically coincide. In the isotropic 
system case, all types of entanglement vanish at the same value of the magnetic field, which is 
small. When we investigated the behavior of the threshold temperatures, we found that they 
exhibit a very similar behavior to the entanglements where the nearest neighbor bipartite 
temperature is much higher than that of the next to nearest neighbor bipartite and multipar- 
tite where the last two coincide asymptotically and the three of them increase monotonically 
with the magnetic field. All the threshold temperatures of the isotropic system vanish at 
the same value of the magnetic field where all the entanglements vanish. It is important 
to emphasis that this means as the temperature increases to certain value, the threshold, 
the multipartite entanglement and the bipartite entanglement of the far from each other 
parts of the system may vanish while the nearest neighbor bipartite entanglements sustain 
up to a much higher temperature. Studying the different systems energy gaps as a function 
of the magnetic field showed that they have great correspondence to the behavior of the 
entanglements and the threshold temperatures, where large characteristic energy gap lead 
to stronger robustness of entanglement and higher threshold temperatures, while vanishing 
energy gap may cause zero threshold temperature, nevertheless the properties of the ground 
state of the system plays a major role in determining the behavior of the entanglement and 
the threshold temperature over the energy gap. This was seen in the isotropic system case, 
which has a ground state that is entangled only below a threshold value of the magnetic 
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field and that caused both the entanglement and the the threshold temperatures to vanish at 
this value regardless of the monotonic increase of energy gap. Furthermore, we have focused 
on examining the persistence of quantum effects at high temperatures where we found that 
the different types of entanglements (specially the bipartite), though would have very small 
values, can be maintained at high temperatures by applying sufficiently high magnetic fields. 
It is interesting in future to investigate the same systems coupled to a dissipative environ- 
ment in presence of thermal excitations to test the interplay of the two environments. Also 
it is important to engineer similar systems with inserted impurities to examine their effect 
to tune the threshold temperatures. Furthermore we would like to investigate the same 
system with larger number of sites to test the system size effect on robustness of thermal 



entanglement and determine threshold temperatures using finite size scaling 
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